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Abstract 

We discuss the results of a global \ 2 analysis of a simple 5*0(10) SUSY GUT 
with _D 3 family symmetry and low energy R parity. The model describes 
fermion mass matrices with 14 parameters and gives excellent fits to 20 
observable masses and mixing angles in both quark and lepton sectors, 
giving 6 predictions. Bi-large neutrino mixing is obtained with hierarchical 
quark and lepton Yukawa matrices; thus avoiding the possibility of large 
lepton flavor violation. The model naturally predicts small 1-3 neutrino 
mixing, with sin #13 ~ 0.05 — 0.06. In this paper we evaluate the predictions 
for the lepton flavor violating processes, /i — > cy, r — > ^7 and r — ► ej 
and also the electric dipole moment of the electron, d e , muon and tau, 
assuming universal squark and slepton masses, mi 6 , and a universal soft 
SUSY breaking A parameter, A , at the GUT scale. We find Br(fi — > e^) is 
naturally below present bounds, but may be observable by MEG. Similarly, 
d e is below present bounds; but is within the range of future experiments. 
We also give predictions for the light Higgs mass (using FeynHiggs). We 
find an upper bound given by < 127 GeV, with an estimated ±3 GeV 
theoretical uncertainty. Finally we present predictions for SUSY particle 
masses in the favored region of parameter space. 



1 Introduction 



In this letter we present results for a global \ 2 analysis of the 50(10) SUSY GUT for 
fermion masses presented in Ref. pp. The model also has a D 3 x [U(l) x Z 2 x Z 3 ] 
family symmetry^] The three families of quarks and leptons are contained in three 
16 dimensional representations of SO (10) {lQ a , I63} with 16 a , a = 1,2 a D3 flavor 
doublet (see Ref. [2] for details on D 3 ). The third family, along with the pair of 
electroweak Higgs doublets, contained in the 10 dimensional representation of 5*0(10), 
are D 3 singlets. Hence only the third generation has a renormalizable Yukawa coupling 
and, as a consequence, we have \ t = = A r = \ Ut Yukawa unification at Mq UT . 
This forces us into the large tan j3 regime and several interesting predictions follow. 
We have derived the consequences of third generation Yukawa unification in several 
papers. In Ref. [5] we demonstrated that in order to fit the low energy values of 
the top, bottom and tau masses (with the typically large, of order 50%, radiative 
corrections to the bottom quark mass) the soft SUSY breaking parameters necessarily 
reside in a very narrow region of the possible parameter space. Hence we have definite 
predictions for SUSY spectra, see [3] and Section H] for more details. In addition, in this 
region of parameter space the light Higgs mass necessarily has a central value of order 
120 GeV. In Ref. [6] we demonstrated that this same minimal 5*0(10) SUSY model 
[MSO10SM] gives the correct abundance of dark matter, fitting the WMAP data, and 
gives observable values for the branching ratio Br(B s — > /i + /i~) with a lower bound 
of order 10" 8 . The dark matter candidate in this model is the LSP, neutralino, which 
predominantly annihilates through a direct s-channel CP odd Higgs, A. In addition, it 
also dominates in the leptonic decay of B s . 

In the present model, all of the above results are retained (with small modifications), 
but in addition we fit the masses and mixing angles of all three families, including neu- 
trino data. The model describes fermion mass matrices with 14 parameters and gives 
excellent fits to 20 observable masses and mixing angles in both quark and lepton sec- 
tors, giving 6 predictions. Both the charged fermion and neutrino mass matrices are 
hierarchical; thus suppressing large flavor violating interactions, even at large tan f3. 
The simple structure of the neutrino sector leads quite naturally to maximal atmo- 
spheric neutrino oscillations and large solar neutrino mixing pQ. We predict a very 
small value for sin^ ~ 0.05 — 0.06. In addition, CP violation in the neutrino sector is 
fixed by the phases in the charged fermion mass matrices. At the same we can easily 
accommodate leptogenesis via non-thermal processes, see for example [TJ. 

2 The Model 

The full superpotential W = W c h,j ermions + W neutr i no for fermion masses and mixing 
angles contains two terms. The first term, resulting in Dirac Yukawa matrices for 
charged fermions and neutrinos, is given by 

.fermions I63 10 I63 + 16 a 10 X a (1) 

+Xa (M x X a + 45 & 16 3 + 45 I 16 a + A 16 a ). 

1 The charged fermion sector of this theory was considered in an earlier paper [5] and the neutrino 
sector of the theory was inspired by the previous analysis by one of us, R.D. [3J. 
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The third family of quarks and leptons is contained in the superfield I63 (transforming 
as a 16 of 50(10)); the first two families are contained in 16 a , a = 1,2 (with explicit 
5*0(10) tranformations and transforming as a D 3 doublet) and the two Higgs doublets 
are contained in 10. The additional fields are an adjoint of 5O(10) (45) and several 
50(10) singlet flavon fields needed to break the full flavor symmetry D 3 x [U(l) x Z 2 x 
Z 3 \. Note, M x = M (l + aX + f3Y) includes 5O(10) breaking vevs in the X and Y 
directions, (fi a , (fi a (D 3 doublets), A (1b singlet) are 5O(10) singlet flavon fields, and 
M, Mo are 50(10) singlet masses. The fields 45, A, <p, <fi are assumed to obtain vevs 
(45) ~ (B — L) Mq (where B, L, Mq is baryon and lepton number and the GUT scale, 
respectively) , A <C M and 

<*>=(£)■<«=(£) (2) 

with 0! > 02- The second term gives large lepton number violating masses for "right- 
handed" neutrinos; necessary for the See-Saw mechanism. We have 

W neutrino = 16(A 2 N a 16 a + A 3 N 3 16 3 ) (3) 
+ I (5 a N a N a + 5 3 N 3 N 3 ) 

where the fields N 3 ; N a a — 1, 2 are 5O(10) singlets and 16 is assumed to break 5O(10) 
to SU (5) via a vev in the right-handed neutrino direction. 

The superpotential, (Eqn. [T]) results in the following Yukawa matrices^] 

e' p 




Y u = \ -e'p ep -e | A 

Y d = 1 -e' e -to I A (4) 
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K = e' uj 31uj \tuj A (5) 



with 



e oc 4> 2 /M; (6) 



e 



oc 0i/M; e' ~ (A/M ) 
1 + a 
1 — 3a ' 



1 + a 

a = — ; p ~ p -C a 



= 20-/(2(7-1). 



2 In our notation, Yukawa matrices couple electroweak doublets on the left to singlets on the right. 
It has been shown in Ref. [8] that excellent fits to charged fermion masses and mixing angles are 
obtained with this Yukawa structure. 
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The Dirac mass matrices are then given by 



v, u (7) 
e, d . (8) 

Consider the neutrino masses. In the three 16s we have three electroweak doublet 
neutrinos (i/ a , v 3 ) and three electroweak singlet anti-neutrinos (u a , u 3 ). In addition, the 
anti-neutrinos get GUT scale masses by mixing with three SO (10) singlets {N a , a = 
1, 2; A^} transforming as a D 3 doublet and singlet respectively. We assume 16 obtains 
a vev, v i6, in the right-handed neutrino direction, and (S a ) = M a for a = 1, 2 (with 
M 2 > M x ) and (S 3 ) = M 3 @ We thus obtain the effective neutrino mass terms given by 

W = um u u + uVN+-NM N N (9) 

with 

/ A 2 \ 

V = v 16 A 2 , M N = diag(M h M 2 , M 3 ). (10) 
\ A 3 / 

The electroweak singlet neutrinos {p, N} have large masses V, ~ Mq. After 
integrating out these heavy neutrinos, we obtain the light neutrino mass matrix given 
by 

M = m v M^ 1 m T u} (11) 
where the effective right-handed neutrino Majorana mass matrix is given by: 

M R = V V T = diag(M Rl , M Ra , M R3 ) , (12) 

with 



rrii = Yi —= sin (3 i 
v2 

mj = Yi—j= cos /3 i 
v2 



Mri = (A 2 v 16 ) 2 /M 2 , M R2 = (A 2 ^le) 2 /^!, Mr 3 = (A 3 u 16 ) 2 /M 3 . (13) 

Defining U e as the 3x3 unitary matrix for left-handed leptons needed to diagonal- 
ize Y e (Eqn. HJ, i.e. F e D = ZTj F e £/? and also U v such that M U u = M D = 
diag(m !/1 , m„ 2 , m„ 3 ), then the neutrino mixing matrix is given by Upmns — U\ U v in 
terms of the flavor eigenstate (v a , a = e, /1, r) and mass eigenstate (z/,, z = 1,2,3) 
basis fields with 

v a = ^2(UpMNs)ai n- (14) 

i 

For Upmns w e use the notation of Ref [9j with 

C12C13 si 2 ci 3 s l3 e~ i5 \ I e iail2 vi 

-s\ 2 c 23 - c 12 s 23 s 13 e lS c i2 c 23 - s 12 s 23 s 13 e lS s 23 c 13 ] ( e lol2/2 u 2 | (15) 
S12S23 - ci 2 c 23 s 13 e lS -c 12 s 23 - s 12 c 23 s 13 e lS c 23 c 13 



3 These are the most general set of vevs for c/> a and S a - The zero vev for <f>i can be enforced with a 
simple superpotential term such as S <j) a <fi a . 
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Finally, we note that this theory is certainly not fundamental with many arbitrary 
symmetry breaking VEVs at the GUT scale. Nevertheless, it has two major features in 
its favor. As a result of the GUT and family symmetries, the Yukawa matrices, which 
are the only observables of the complicated GUT physics, have fewer parameters than 
low energy observables. Hence this theory is predictive. Secondly, the model has the 
advantage that it self-consistently fits the low energy data, and thus, at the very least, 
it is an excellent phenomenological ansatz for fermion masses. Thus it can be tested 
via additional low energy flavor violating processes. 

3 Global x 2 Analysis 

Yukawa matrices in this model are described by seven real parameters {A, e, e, a, 
p, e', £} and, in general, four phases $g, Therefore, in the charged 

fermion sector we have 11 parameters to explain 9 masses and three mixing angles 
and one CP violating phase in the CKM matrix, leaving us with 2 predictionsjj Note, 
these parameters also determine the neutrino Yukawa matrix. Finally, our minimal 
ansatz for the right-handed neutrino mass matrix is given in terms of three additional 
real parameters^, i.e. the three right-handed neutrino masses. At this point the three 
light neutrino masses and the neutrino mixing matrix, Upmnsi (3+4 observables) are 
completely specified. Altogether, the model describes 20 observables in the quark and 
lepton sectors with 14 parameters, effectively having 6 predictions^ 

In addition to the parameters describing the fermion mass matrices, we have to input 
three parameters specifying the three gauge couplings: the GUT scale M G defined as 
the scale at which a,\ and a 2 unify, the gauge coupling at the GUT scale a G , and the 
correction e 3 to a% {Mq) necessary to fit the low energy value of the strong coupling 
constant. Finally we have to input 7 supersymmetry parameters given by - M1/2, 
a universal gaugino mass; mio, a universal Higgs mass; rriiQ, a universal squark and 
slepton mass; A , a universal trilinear coupling; a small Higgs mass splitting parameter, 
Amu = l/2(m 2 Hd — Tn 2 H ^)/m\ Q ] the supersymmetric Higgs mass parameter \i{Mz) and 
the ratio of the two Higgs VEVs, tan (3. 

We have also imposed some physically motivated constraints on the x 2 analysis. We 
demand a lower bound on the lightest stop mass given by m 4 ~ = 500 GeV. Lower values 
of mi actually give even better fits. On the other hand, a chargino-stop loop gives the 
dominant SUSY contribution to the process b — > S7 and lighter stop values make it 
difficult to fit this process. In addition we fix the CP odd Higgs mass, rriA = 700 GeV. 
Lower values would result in a branching ratio Br(B s — > [i + which approaches the 
experimental lower bound, see Ref. [6]. Our results however are not sensitive to this 
latter constraint. Further discussion of the former constraint is given in Sect. HI 

4 Of course, in any supersymmetric theory there is one additional parameter in the fermion mass 
matrices, i.e. tan/3. Including this parameter, there is one less prediction for fermion masses, but then 
(once SUSY is discovered) we have one more prediction. This is why we have not included it explicitly 
in the preceding discussion. 

5 In principle, these parameters can be complex. We will nevertheless assume that they are real; 
hence there are no additional CP violating phases in the neutrino sector. 

6 Note, the two Majorana phases are in principle observable, for example, in neutrinoless double- 
beta decay [TU], however, the measurement would be very difficult (perhaps too difficult [H]). If 
observable, this would increase the number of predictions to 8. 
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All the parameters (except for fi(Mz)) are run from the GUT scale to the weak 
scale (Mz) using two (one) loop RGEs for dimensionless (dimensionful) parameters. 
At the weak scale, the SUSY partners are integrated out leaving the two Higgs doublet 
model as an effective theory. We require proper electroweak symmetry breaking. More- 
over, the full set of one loop, electroweak and SUSY, threshold corrections to fermion 
mass matrices, as well as to the three precision electroweak observables (G M , a^M an d 
a s (Mz)) are calculated at MzQ Below Mz we use 3 loop QCD and 1 loop QED RG 
equations to calculate light fermion masses and oiem- More details about the analysis 
can be found in [5] or j(i El In addition, we self-consistently include the contributions 
of the right-handed neutrinos to the RG running between the GUT scale and the mass 
of the heaviest right-handed neutrino [12]. 

The x 2 function is constructed from observables given in Table [2j We have used 
the top quark mass from the PDG reviews of particle properties 2005 update by T. 
M. Liss and A. Quadt@ Note that we have included several redundant observables in 
the quark sector. We do this because quark masses are not known with high accuracy 
and different combinations of quark masses usually have independent experimental and 
theoretical uncertainties. Thus we include three observables for the charm and bottom 
quark masses: the MS running masses (m c (m c ), m&(mb)) and the difference in pole 
masses M& — M c obtained from heavy quark effective theory. The same is true for 
observables in the CKM matrix. For example, we include Vtd and the two CP violating 
observables ex and the value for sin (2/3) given by the world average measured via the 
process B —>■ J ftp K s |13j . However we have doubled the error to take into account 
the significant difference between the BaBar and Belle central values. We thus have 16 
observables in the quark and charged lepton sectors. We use the central experimental 
values and one sigma error bars from the particle data group [9|. However, in the case 
that the experimental error is less than 0.1% we use a = 0.1% due to the numerical 
precision of our calculation. 

At present only four observables in the neutrino sector are measured. These are 
the two neutrino mass squared differences, Am^ and Am^, and two mixing angles, 
sin 2 9 12 and sin 2 # 2 3- For these observables we use the central values and 2a errors from 
Ref. [H]. The other observables: neutrino masses, 1-3 mixing angle and the phase of 
the lepton mixing matrix are predictions of the model. In addition, the new feature 
of this paper is the predictions for several lepton flavor violating processes and lepton 
electric dipole moments. 



4 Results 

Let us now discuss our results. We performed the global \ 2 analysis for values of the 
soft SUSY breaking scalar mass at Mq given by m 16 = 3, 4 and 5 TeV. However, we 
only present the results for the latter two cases. Good fits prefer the region of SUSY 

In a top down analysis, the DR value of sin 2 0w{Mz) is obtained directly via RG running <xn t 2) 
from the GUT scale. Then with the calculated value of Mz and all one loop threshold corrections 
included in Ar, we obtain the observed value of G M [4]. 

8 The only difference is that in the present analysis we include all three families of fermions. 

9 This value of the top quark mass (172.7 ± 2.9 Gev/c 2 ) is not so different from the most recent 
value (172.5 ± 1.3 stat ± 1.9 syst GeV /c 2 ) from CDF II and DZero at the Tevatron [15]. 
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parameter space characterized b; 



H, M1/2 < m 16 (16) 
-A ~ v^mio ~ 2mi 6 . 

This is required in order to fit the top, bottom and tau masses when the third gen- 
eration Yukawa couplings unify [5j. Note, the three input parameters (/i, M 1 / 2 , m 16 ) 
are not varied when minimizing y 2 . As a consequence of the relations, (Eqn. [IB]) , we 
expect heavy first and second generation squarks and sleptons, while the third genera- 
tion scalars are significantly lighter (with the stop generically the lightest). In addition, 
charginos and neutralinos are typically the lightest superpartners. We predict values of 
tan/3 ~ 50 and a light Higgs with mass, rrih < 127 GeV (with a theoretical uncertainty 
±3 GeV). The specific relations between the SUSY breaking parameters also lead to 
an interesting prediction for the process B„ — » with branching ratio in the region 
currently being explored at the TevatronJUj Furthermore, the neutralino relic density 
obtained for our best fit parameters is consistent with WMAP data [6J and direct neu- 
tralino detection is possible in near future experiments. Finally, this region maximally 
suppresses the dimension five contribution to proton decay [16J and suppresses SUSY 
flavor and CP violation in general. For more information on the SUSY and Higgs spec- 
tra and related phenomenology in this region of SUSY breaking parameter space, see 
Refs. [5] and [6J. 

In Figs. [1] and [2] we present contours of constant y 2 f° r m w — 4 and 5 TeV with 
rriA = 700 GeV, as a function of /i, My 2- The best fits are obtained for small values of 
M\/2 < 300 GeV (where the lower bound on My% is determined by the experimental 
bound on the chargino mass, m x + > 104 GeV). We find that the value of x 2 decreases 
as mi6 increases. This is solely due to the lower bound of 500 GeV on the lightest stop 
mass and the resultant difficulty in fitting the bottom quark mass with heavier stop. 

At this point a brief aside is necessary. In our analysis, we have not evaluated 
several significant pieces of data. These include the branching ratios Br(b — * 37), 
BriB — > X s l + / _ ), and B s — B s mixing. These also provide significant constraints on 
the theory. However, we have used the code of T. Blazek [TTJ, [S] to check our analysis 
and also evaluate the branching ratio Br(b — > 57). This process is enhanced at large 
tan (3. The dominant SUSY contribution comes from the chargino-stop loop. We find 
that the most significant constraint is a lower bound on the lightest stop mass of order 
500 GeV. We have thus imposed this bound on the stop mass by introducing a penalty 
to x 2 ■ The best fit for the branching ratio Brib — > 37) is then fit with the minimal stop 
mass. Note, we find that x 2 increases as the lower bound on the stop mass increases. 
This is due to the fact that a good fit for m& prefers a light stop [5j [6J . In addition, 
as the lower bound on the stop mass increases, we find it necessary to increase m 16 . 
For example, with the light stop mass, mi = 300 GeV, we find good fits with m 16 = 3 
TeV [5J. Now with m ; - = 500 GeV, good fits, with x 2 < 8, are only obtained with 

10 In addition, the best fit requires a non-universal Higgs masses at the GUT scale with Am# = 
l/2(rn? H — rn 2 Hu )/rn\ ~ .07. Note this is significantly smaller than needed in the past [5]. That is 
because the RGE running of neutrino Yukawas from Mq to the heaviest right handed neutrino has 
been included self-consistently. As noted in [S], such running was a possible source for Higgs splitting. 
Evidently it can not be the only source. 

11 This process is sensitive to the CP odd Higgs mass, itia, which can be adjusted in theories with 
non-universal Higgs masses. 
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Figure 1: Contours of constant x 2 f° r m w = 4 TeV and itia = 700 GeV. The yellow (very 
light) shaded region at the bottom and top (and the region bounded by the extended solid 
boundary line) has x 2 > 8. The blue (light shaded) region on the left (and below the extended 
dotted line) is excluded by m x + < 104 GeV and the green (darker shaded) region is excluded 
by the Higgs mass bound mh < 111 GeV. 



7^16 > 4 TeV. As a final note, for a light stop (m 4 ~ = 500 GeV) the Wilson coefficient of 
the dominant operator for the process b — > 57, O7, has the opposite sign in the MSSM 
than for the standard model, i.e. Cj(MSSM) ~ —Cj(SM) |18j . Recent measurements 
of the branching ratio Br(B — >• X s l + l~) [19] suggest that the same sign is preferred, 
Cj(MSSM) ~ Cr(SM); while the latest data on the forward-backward asymmetry 
does not seem to distinguish these two possibilities [20]. Clearly a more detailed \ 2 
analysis including all these processes would be necessary to better test the theory. This 
is however not the focus of the present paperl^l 

In Figs. [3] and H] we present contours of constant light Higgs mass for the case 
rriiQ = 4 and 5 TeV. There is not much difference in the range of light Higgs mass in 
the two cases. We find an upper bound on the Higgs mass given by < 127 GeV. 
In our analysis, we use the output of our RG running as input to FeynHiggs [21] to 
obtain the Higgs pole mass at two loops. Note, in the region of parameter space with 
\A \ > m 16 the radiative corrections to the Higgs mass are significant, i.e. the two 
loop correction, using FeynHiggs, is of order 30 GeV. One might then worry that the 
theoretical uncertainty in the light Higgs mass is just as big. However, Heinemeyer 
[22] (see sec 2.5), estimates the uncertainties in the light Higgs mass from yet-to-be- 
calculated two-loop corrections and higher to be at most 3 GeV. We have thus taken 

12 We have also not calculated (<? — 2) M in this paper. However in a previous analysis [5] we found 
typical values of (g — 2) M < 3 x 10~ 10 . 
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Figure 2: Contours of constant x 2 f° r m i6 = 5 TeV and tua = 700 GeV. The yellow (very 
light) shaded region at the bottom and top (and the region bounded by the extended boundary 
line) has x 2 > 6. The blue (light shaded) region on the left (and below the extended dotted 
line) is excluded by rn x + < 104 GeV and the green (darker shaded) region is excluded by the 
Higgs mass bound mh < 111 GeV. 



±3 GeV as the estimated total theoretical uncertainty in the light Higgs massP*! 

4.1 Lepton flavor violation and electric dipole moments 

Let us now focus on our results for lepton flavor violating (LFV) processes lj — ► la 
and charged lepton electric dipole moments (EDMs) in this theory. We start with 
universal squark and slepton masses and a universal A parameter at the GUT scale. 
This is thus a GUT version of minimal flavor violating boundary conditions, giving 
minimal flavor violation at low energies. Thus the dominant contribution to lepton 
flavor violation results from the RG running of slepton masses and the effect of neutrino 
Yukawa couplings on this running from the GUT scale to the heaviest right-handed 
neutrino Majorana mass of order 10 14 GeV. See for example the seminal paper on this 
subject |24j . There is also ample literature regarding LFV and EDMs, see for instance 

13 Our results for the light Higgs mass are about 5 GeV higher than in our previous analyses [SI [6]. 
This is due to the fact that we now use FeynHiggs to calculate the Higgs mass, where as in previous 
papers we used an effective potential analysis. The disagreement of the effective potential method 
with FeynHiggs (a perturbative analysis) is well known, see for example 23J. In addition one sees that 
the light Higgs mass decreases as M1/2 increases. This fact is completely due to the global y 2 analysis 
and the need to fit the bottom quark mass starting with 3rd family Yukawa unification. For a detailed 
discussion of this effect, see the second reference in [BJ. 
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Figure 3: Contours of constant light Higgs mass nih for ra\Q = 4 TeV and rriA = 700 GeV. 
The shaded regions are the same as in Fig. 1. 



[251 126] for LFV and (23, [28] for EDMs. Therefore we simply quote the results of [25] 
and [27] here and refer the reader to those references for more detail. 

Following the notation of [25], the effective Lagrangian £ relevant for the decay 
lj -> la is 

£ = -^m h u l( j a p(ALp L + A R P R ) Uj F a(3 } 

where e is the electric charge, m\ j is the mass of the decaying lepton, Pr/l is the chirality 

projection operator, Ui and Uj are Dirac spinors describing U and lj, respectively. A^ R 
is obtained by calculating Feynman diagrams depicted in Figure [5] at one loop, and 
found in [2"5] . 

The decay amplitude is given by 

T = iee^mt.Uiip - q)a a ^(A^P L + A R P R ) Uj (p), 

where e* is the photon polarization vector. Then, the decay rate is 

T(l J ^ln) = ^ml(\A^ + \A R \ 2 ). 

On the other hand, the lepton electric dipole moment is defined as the coefficient 
of the effective Lagrangian £ of the form 

£ = -^d h u i a at3 'y 5 u i F al3 . 
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Figure 4: Contours of constant light Higgs mass nih for ra\Q = 5 TeV and rriA = 700 GeV. 
The shaded regions are the same as in Fig. 2. 




Figure 5: One-loop Feynman diagrams relevant for the decay of lj — ► ^7. (a) involves 
charginos x° A and sneutrinos u a , and (b) involves neutralinos and sleptons l a in the loop 
with A = 1, 2 for charginos, A = 1, 2, 3, 4 for neutralinos, a = 1, 2, 3 and a = 1, . . . , 6. 



Let us write 

^ = df + df , 

where df h and <i"* are contributions to the EDM from loops in Figure [5] (a) and (b) 
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with replacing L by k, respectively. Then we find 



2 3 

jch _ e t™ (r L ® r R ( l > \ 3-4x Aa + x 2 4a + 21nx Aa 

^ --T^A,l. Im ^ C ^ ^ 2(l-x Aa )3 > 

^ _ e A A L(0 o (0 m^o l-x Aa + 2x A alnx Aa 

Here Cfjf^ , C!^P , and are read from vertices shown in Figure [6] and the 

expression for them is given in [25]. m^,+ is chargino mass, mp a is sneutrino mass, m^o 
is neutralino mass, is selectron mass, x Aa = m 2 + /m? a and x Aa = m| /m~ . 



Xa 



Xa 



^^P R + C^P L ) 



(b) 



Figure 6: Vertices (a) define c£j£ and C^?? and vertices (b) iVj£ and iV^ with A = 1,2 
for charginos, j4 = 1, 2, 3, 4 for neutralinos, and a = 1, . . . , 6. 



Table [T] shows the current limits on various lepton flavor violating processes and 
EDMs with an estimated sensitivity for future experiments. See for example [42J for a 
summary of the current and future experimental status on searches for LFV and the 
muon EDM. 





Current limit 


Expected accuracy 


Br(n — > ej) 
Br(r — > ey) 
Br(r — > fi'-f) 
Br(n — ► 3e) 
5r(r -> 3/) 
Br(fiTi — > eTi) 


< 1.2 x 10" 11 [29] 

< 1.1 x 10" 7 [3T] 

< 6.8 x 10" 8 [32] 

< 1.0 x 10~ 12 [33] 

< 1 - 3 x 10~ 7 [31] 

< 1.7 x 10~ 12 [35] 


5 x 10" 14 [30] 
10" 18 [36] 


c? e [e • cm] 
[e ■ cm] 
d r [e • cm] 


< 1.6 x 10" 27 [37] 

< 10" 18 [39] 

< 10" 16 [H] 


10" 31 [38] 

io~ 24 - to- 25 do] 



Table 1: Current upper limit on various LFV processes and lepton EDMs. The limits 
on the LFV processes are all at 90 % CL. Estimated accuracy for the future experiment 
is given for // — > ej, fie conversion, electron EDM and muon EDM. 

We note that the MEG collaboration [30] (searching for \i — > c-f with a sensitivity of 
order > 5 x 10~ 14 ) should start taking data by September 2006 and may have significant 
results by 2008. It will be an excellent test for any new physics beyond the standard 
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model; and in particular, the DR model pp. Let us also note here that we can calculate 
Br(U 31 j) to a good approximation by using 

Br(l t -> -3n{ g ml 4 J ' 

which has been verified in Ref. [13] for all values of tan (3. In particular 

Br(r -> 3/x) ^ 1 fir(r -> 3e) _ 1 Br(/x -> 3e) _ 1 
5r(r -»• ~ 440' Br(r -» i^) ~ 94' 5r(/i -» i^) ~ 162' 

This means that if we satisfy the constraint from Brili — > ^7), we automatically also 
satisfy the constraint from -Br(£j — > 3Zj- ) . 

In Fig. [7] we plot contours of constant branching ratio Br (fx — > ej) for mig = 4 
TeV. The prediction is significantly below the present experimental bounds. Moreover, 
comparing this result with the future sensitivity of the MEG experiment [30] (Br(fi — ► 
ej) > 5 x 10~ 14 ) we find that our prediction is below the MEG sensitivity in most of 
the parameter space. Note, however, the narrow region in the upper right hand corner 
with x 2 < 8 which is within the sensitivity of the MEG experiment. In Fig. [S] we 
present results for mie = 5 TeV. Of course, larger scalar masses suppress the branching 
ratio, so that now the entire allowed region is below the projected MEG sensitivity. 
The results for the decays r — > cy and r — > /ry are given in Fig. [91 Unfortunately 
the results are significantly below the present bounds and we are not aware of any 
experiments to significantly improve these bounds. 



Br(n->e, y)x10 14 




100 150 200 250 300 350 400 450 500 

M 1/2 [GeV] 

Figure 7: Contours of constant branching ratio Br([i — > ej) x 10 14 for m\Q = 4 TeV and 
niA = 700 GeV. The shaded regions are the same as in Fig. 1. Note the narrow region in the 
upper right hand corner with x 2 < 8 which is within the sensitivity of the MEG experiment. 
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Br(ji^e, y)x10 15 




100 150 200 250 300 350 400 450 500 

M 1/2 [GeV] 

Figure 8: Contours of constant branching ratio Br(fi — > erf) x 10 15 for miQ = 5 TeV and 
rriA = 700 GeV. The shaded regions are the same as in Fig. 2. Note all points are below the 
sensitivity of the MEG experiment. 



We have also evaluated the predictions for the electric dipole moment of the electron, 
muon and tau. In Figs. [10] and [11] we present the results for the electric dipole moment 
of the electron. Note, in both cases, the entire region is below the present bounds, and 
also within the projected sensitivity of future experiments [38]. In Fig. [121 we present 
the results for the electric dipole moments of the muon and tau. In all cases the results 
are below the present bounds and for the result is below the projected sensitivity of 
future experiments [4"0] . 

In Figs. [TBI and [T4l we evaluate the neutrino mixing angle sin 2 #13. Recall that 
measuring this mixing angle is the goal of several future reactor and long baseline 
neutrino experiments. Moreover, a sufficiently large value for sin 2 # 13 is needed in order 
to have the possibility of observing CP violation in neutrino oscillations. The value of 
sin 2 #i3 is somewhat sensitive to the value of mie; with a central value changing from 
sin 2 #13 ~ 0.0030±0.0007 for m 16 = 4 TeV to sin 2 6 13 ~ 0.0024±0.0004 for m ie = 5 TeV 
(where the uncertainty corresponds to varying over the range for fi, M\i2 with x 2 <9> 
or < 6 in the two cases). Note, while sin 2 # 13 is relatively insensitive to varying /1, M\ii\ 
the CP violating parameter sin 5, on the other hand, is quite sensitive. We find that 
sin 5 can vary between 0.1 and 1.0 for different values of /x, M1/2, and as a result the CP 
violating Jarlskog parameter J ranges between 0.0013 and 0.013. CP violation in the 
latter case may be observable at long baseline experiments. For example, the JPARC- 
SK experiment has a potential sensitivity to sin 2 26 Vi < 1.5 x 10 -3 and 5 ~ ±20° and 
a comparable sensitivity is expected from the "Off-axis NUMI" proposal [44J . 

In Tables [2] and [3] we present the x 2 fit for a particular point in SUSY parameter 
space for mi 6 = 4 and 5 TeV, respectively. The points give a value of x 2 — 7.65 and 
4.99. The former is acceptable while the latter is quite good. In the table caption we 
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Figure 9: Contours of constant branching ratio Br(r — ► e~f) x 10 12 (upper) and Br(r — > 
^7) x 10 11 (lower) for mi6 = 4 TeV and toa = 700 GeV. The shaded regions are the same as 
in Fig. L 



present the input data at the GUT scale. We also show the heavy Majorana neutrino 
masses (roughly 10 10 , 10 12 , 10 14 GeV) responsible for the See-Saw mechanism and the 
light neutrino masses. 

Note, that the pull from and Mt> — M c is significantly lower for m 16 = 5 TeV than 
for mi6 = 4 TeV. This accentuates the "tug of war" between the gluino and chargino 
loop contributions to the bottom quark mass at large tan f3. The light quark mass 
ratio rrid/m s is difficult to fit and contributes significantly to the pull in both cases. 
This mass ratio is particularly sensitive to the Georgi-Jarlskog ansatz relating first and 
second generation quark and lepton masses. The conflict here is with the very low 
value of the strange quark mass, of order 105 MeV, preferred by lattice gauge theory 
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d e [ecm] x 10' 




100 150 200 250 300 350 400 450 500 

M 1/2 [GeV] 

Figure 10: Contours of constant electric dipole moment of the electron, d e [e ■ cm] x 10 29 , for 
m\Q = 4 TeV and rriA = 700 GeV. The shaded regions are the same as in Fig. 1. Note the 
entire region is within the sensitivity of future experiments. 



calculations. Finally, we note that in a previous analysis [I] sin 2/3 contributed a value 
of 1.5 to the pull. However the recent Belle data gives a significantly smaller central 
value for sin 2(3 and now the fit is significantly improved. 

We present the additional predictions for squark, slepton and Higgs masses, at these 
two points in SUSY parameter space, in Table HJ and for neutrino masses and mixing 
parameters, in Table El We have given the value for the effective mass parameter 
observable in neutrinoless double beta decay 



y.u 2 

Z_/« e 



E 



U, 



(17) 



(where a[ = ai + 25, i = 1, 2 [32]). It is predicted to be of order 2 x 10 -4 eV which is 
too low to see in near- future experiments [3HJ |33]- We also give the effective electron- 
neutrino mass observable, relevant for the analysis of the low energy beta decay of 
tritium. This mass parameter is unaffected by Majorana phases and is predicted to be 
an order of magnitude larger. The observable, 



m 



eff 



E 



1/2 



m, 



is predicted to be of order 6 x 10~ 3 eV. The current experimental limit is vn e JJ < 2.5 
eV with the possibility of future experiments, such as KATRIN, reaching bounds on 
the order of 0.35 eV [33]. Unfortunately, both mass parameters may be unobservable 
by presently proposed experiments. Finally, in Table El we present the predictions for 
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d e [e cm] x 1 0' 




100 150 200 250 300 350 400 450 500 

M 1/2 [GeV] 

Figure 11: Contours of constant electric dipole moment of the electron, d e [e ■ cm] x 10 29 , for 
mi6 = 5 TeV and tha = 700 GeV. The shaded regions are the same as in Fig. 2. Note the 
entire region is within the sensitivity of future experiments. 



lepton flavor violation and the electric dipole moments at the same points in SUSY 
parameter space. 

5 Summary and Conclusions 

In this paper we have performed a global x 2 analysis on a well-motivated, phenomeno- 
logically acceptable minimal SO(10) SUSY GUT with a D% family symmetry. The 
most stringent constraint comes from assuming Yukawa coupling unification for the 
third family of quarks and leptons. The \ 2 contours as functions of fi and Mi/ 2 for 
m 16 = 4 and 5 TeV are given in Figs. 1 and 2. We find acceptable solutions with 
X 2 < 8 (6) for miQ = 4 (5) TeV, respectively. We find the light Higgs mass, found using 
FeynHiggs, has an upper bound given by m h < 127 GeV (see Figs. [3]and|I|). The 
additional predictions for the SUSY spectrum, and neutrino masses and mixing angles 
at two particular points in SUSY parameter space are given in Tables |4] and [5j 

In addition to the global \ 2 analysis, we focused on obtaining the rates for several 
lepton flavor violating processes and also for the charged lepton electric dipole moments. 
We calculated the branching ratios for the lepton flavor violating processes Br(fi — > e^) 
(Figs. [7| and [8]) and r — > ey and r — > /ry (Fig. [9]). There is only a narrow region in 
the upper right hand corner (Fig. [7]) with \ 2 ^ 8 which is within the sensitivity of the 
MEG experiment. We have also evaluated the electric dipole moments of the electron 
(Figs. dn]and[TT]) and muon and tau (Fig. [T2]). In all cases the results are below the 
present experimental bounds, and for the result is below the projected sensitivity 
of future experiments [40]. However, in both cases, the entire region is within the 
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Table 2: The fit for for fermion masses and mixing angles at one particular 
point in SUSY parameter space defined by mi 6 = 4 TeV, /i = 300 GeV and 
Mi/a = 200 GeV. 

Initial parameters: (1/«g, M g , e 3 ) = (24.84, 3.30 x 10 16 GeV, -3.57 %), 

(A, Ae, cr, Ae, p, Ae', Ae£) = (0.62, 0.030, 0.87, 0.0063, -0.059, -0.0021, 0.0040), 

($ a , $ e ~, $ p , $ 5 ) = (0.637, 0.453, 0.709, 3.609) rad, 

(mi 6 , Mi/ 2 , A), M M z)) = (4000, 200, -7809.1, 300) GeV, 

{{m H Jm m )\ (m Hu /m ie ) 2 , tan/5) = (1.91, 1.61, 50.34) 

(M R3 , M R2 , M Rl ) = (4.6 x 10 13 GeV, -8.1 x 10 11 GeV, 1.1 x 10 10 GeV) 



Observable 


Data (a) 


1 neory 


T ) 11 

rull 


(masses in GeV) 








x 10 & 


1.16637 (0.1%) 


1.16638 


< 0.01 


—i 

a EM 


137.036 (0.1%) 


137.035 


< 0.01 


a s (M z ) 


0.1187 (0.002) 


0.1174 


n o t 

0.37 


M t 


1 TO T / O r\ \ 

172.7 (2.9) 


1 TO 1 1 

173.11 


r\ AO 

0.02 


m b {M b ) 


4. ZD ^U.ZOJ 


A AO 

4.4y 


U.y4 


M b — M c 


3.4 (0.2) 


3.61 


1.16 


m c (m c ) 


1.2 (0.2) 


1.16 


0.03 


m s 


0.105 (0.025) 


0.107 


0.01 


m d /m s 
Q- 2 x 10 3 


0.0521 (0.0067) 


0.0638 


3.09 


1.934 (0.334) 


1.815 


0.12 


M T 


1.777 (0.1%) 


1.777 


< 0.01 


M M 


0.10566 (0.1%) 


0.10566 


< 0.01 


M e X 10 3 


0.511 (0.1%) 


0.511 


< 0.01 


V us 


0.22 (0.0026) 


0.2193 


0.06 


v cb 


0.0413 (0.0015) 


0.0410 


0.03 


v ub 


0.00367 (0.00047) 


0.00316 


1.15 


V td 


0.0082 (0.00082) 


0.00824 


< 0.01 




0.00228 (0.000228) 


0.00234 


0.08 


sin(2/3) 


0.687 (0.064) 


0.6435 


0.46 


Am^ x 10 3 


2.3 (0.6) 


2.382 


0.01 


Amli x 10 5 


7.9 (0.6) 


7.880 


< 0.01 


sin 2 6*i2 


0.295 (0.045) 


0.289 


0.01 


sin 2 #23 


0.51 (0.13) 


0.532 


0.03 


TOTAL x 2 






7.65 
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Table 3: The fit for for fermion masses and mixing angles at one particular 
point in SUSY parameter space defined by mi 6 = 5 TeV, /i = 300 GeV and 
Mi/a = 200 GeV. 

Initial parameters: (1/«g, M g , e 3 ) = (24.90, 3.29 x 10 16 GeV, -3.45 %), 

(A, Ae, cr, Ae, p, Ae', Ae£) = (0.63, 0.030, 0.77, 0.0070, -0.054, -0.0022, 0.0035), 

($ a , $ e ~, $ p , $ 5 ) = (0.643, 0.410, 0.692, 3.618) rad, 

(mi 6 , Mi/ 2 , A), »(M Z )) = (5000, 200, -9918.9, 300) GeV, 

{{m H Jm m )\ (m Hu /m ie ) 2 , tan/5) = (1.91, 1.62, 50.53) 

(M R3 , M R2 , M Rl ) = (6.1 x 10 13 GeV, -8.6 x 10 11 GeV, 9.6 x 10 9 GeV) 



Observable 


Data (a) 


1 neory 


T ) 11 

rull 


(masses in GeV) 








x 10 & 


1.16637 (0.1%) 


1.16638 


< 0.01 


—i 

a EM 


137.036 (0.1%) 


137.036 


< 0.01 


a s (M z ) 


0.1187 (0.002) 


0.1178 


0.18 


M t 


172.7 (2.9) 


1 TO /"* /I 

173.64 


0.10 


m b {M b ) 


/I OK IC\ OK\ 
4. ZD ^U.ZOJ 


/I on 


U.Uo 


M b — M c 


3.4 (0.2) 


3.49 


0.23 


m c (m c ) 


1.2 (0.2) 


1.08 


0.34 


m s 


0.105 (0.025) 


0.113 


0.12 


m d /m s 
Q- 2 x 10 3 


0.0521 (0.0067) 


0.0629 


2.60 


1.934 (0.334) 


1.824 


0.10 


M T 


1.777 (0.1%) 


1.777 


< 0.01 


M M 


0.10566 (0.1%) 


0.10566 


< 0.01 


M e X 10 3 


0.511 (0.1%) 


0.511 


< 0.01 


V us 


0.22 (0.0026) 


0.2195 


0.03 


v cb 


0.0413 (0.0015) 


0.0413 


< 0.01 


v ub 


0.00367 (0.00047) 


0.00325 


0.77 


V td 


0.0082 (0.00082) 


0.00812 


< 0.01 




0.00228 (0.000228) 


0.00231 


< 0.01 


sin(2/3) 


0.687 (0.064) 


0.6511 


0.31 


Am^ x 10 3 


2.3 (0.6) 


2.408 


0.03 


Amli x 10 5 


7.9 (0.6) 


7.880 


< 0.01 


sin 2 6*i2 


0.295 (0.045) 


0.288 


0.01 


sin 2 #23 


0.51 (0.13) 


0.537 


0.04 


TOTAL x 2 






4.99 
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Table 4: Predictions for SUSY and Higgs spectra for the fit given in Tables [2] and [3] in 
units of GeV. 



Particle 


rriiQ = o lev 


rriiQ = 4 lev 


n 


120 


124 


ti 


oyy 


oyy 


AO 


( uu 


oyy 


H + 


I uu 


7D1 

1 U X 


X? 


81 


80 


xl 


151 


150 


xt 


151 


150 


9 


605 


597 


ti 


498 


496 


h 


902 


846 


Tl 


1686 


1421 



Table 5: Predictions for neutrino masses, sin 2 #13 and CP violation in the lepton sector 
for the fit given in Tables [2] and [31 





mi6 = 5 TeV 


mi6 = 4 TeV 


m V3 (eV) 


0.0492 


0.0489 


m U2 (eV) 


0.0097 


0.0096 


m Ul (eV) 


0.0041 


0.0036 


sin 2 6»i3 


0.0025 


0.0037 


J 


0.0013 


0.013 


sin 5 


0.119 


0.996 


ai (rad) 


-2.974 


-1.771 


a 2 (rad) 


0.136 


1.315 


(mpp) (eV) 


0.00010 


0.00038 


rrfJJ (eV) 


0.0067 


0.0067 


ei 


0.92 x 10~ 7 


1.61 x 10~ 7 



Table 6: Predictions for branching ratios for lepton flavor violating processes and the 
electric dipole moment of leptons for the fit given in Tables [2] and [31 





m 16 = 5 TeV 


m 16 = 4 TeV 


£>r(/i — > cy) 
Br{j — > ej) 
Br{r — > /i7) 
d e [e ■ cm] 
[e ■ cm] 
d T [e ■ cm] 


4.69 x HT 15 
1.26 x 10~ 12 
6.13 x 10~ u 
9.15 x 10~ 30 
-3.58 x 10~ 28 
-1.54 x 10~ 27 


1.40 x 10" 14 
2.40 x 10~ 12 
1.22 x 10~ 10 
2.43 x 10~ 29 
-7.93 x 10^ 28 
-2.64 x 10" 27 
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Figure 12: Contours of constant electric dipole moment of the muon, d^e-cm] x 10 28 , (upper) 
and the tau, d T [e ■ cm] x 10 27 , (lower) for m\§ = 4 TeV and vtia = 700 GeV. The shaded 
regions are the same as in Fig. 1. Note, the entire region is below present bounds and for d^ 
the result is below the projected sensitivity of future experiments. 



projected sensitivity of future experiments for d e [38J. The results for the lepton flavor 
violating processes and electric dipole moments at the same two particular points in 
SUSY parameter space are given in Table [6j 
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